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Abstract 

We generalize the dynamical-mean field (DMFT) approximation by including into the DMFT 
equations some length scale via a (momentum dependent) "external" self-energy S^. This external 
self-energy describes non-local dynamical correlations induced by short-ranged collective SDW- 
like antiferromagnetic spin (or CDW-like charge) fluctuations. At high enough temperatures these 
fluctuations can be viewed as a quenched Gaussian random field with finite correlation length. 
This generalized DMFT+Sk approach is used for the numerical solution of the weakly doped one- 
band Hubbard model with repulsive Coulomb interaction on a square lattice with nearest and next 
nearest neighbour hopping. The effective single impurity problem in this generalized DMFT+Sk 
is solved by numerical renormalization group (NRG). Both types of strongly correlated metals, 
namely (i) doped Mott insulator and (ii) the case of bandwidth W < U (U — value of local 
Coulomb interaction) are considered. Densities of states, spectral functions and ARPES spectra 
calculated within DMFT+Sk show a pseudogap formation near the Fermi level of the quasiparticle 
band. 

PACS numbers: 71.10.Fd, 71.10.Hf, 71.27+a, 71.30.+h, 74.72.-h 



I. INTRODUCTION 



Among the numerous anomalies of the normal phase of high-temperature superconduc- 
tors the observation of a pseudogap in the electronic spectrum of underdoped copper oxides^ 
is especially interesting. Despite continuing discussions on the nature of the pseudogap, the 
preferable "scenario" for its formation is most likely based on the model of strong scattering 
of the charge carriers by short-ranged antiferromagnetic (AFM, SDW) spin fluctuations 2 ^. 
In momentum representation this scattering transfers momenta of the order of Q = (-, -) 
(a — lattice constant of two dimensional lattice). This leads to the formation of struc- 
tures in the one-particle spectrum, which are precursors of the changes in the spectra due 
to long-range AFM order (period doubling). As a result we obtain non-Fermi liquid like 
behavior (dielectrization) of the spectral density in the vicinity of the so called "hot spots" 
on the Fermi surface, appearing at intersections of the Fermi surface with antiferromagnetic 
Brillouin zone boundary 2 . 

Within this spin-fluctuation scenario a simplified model of the pseudogap state was 
studied 2 ^*** 1 under the assumption that the scattering by dynamic spin fluctuations can be 
reduced for high enough temperatures to a static Gaussian random field (quenched disor- 
der) of pseudogap fluctuations. These fluctuations are defined by a characteristic scattering 
vector from the vicinity of Q, with a width determined by the inverse correlation length of 
short-range order k = 

Undoped cuprates are antiferromagnetic Mott insulators with U 3> W (U — value of 
local Coulomb interaction, W — bandwidth of non-interacting band), so that correlation 
effects are very important. It is thus clear that the electronic properties of underdoped 
(and probably also optimally doped) cuprates are governed by strong electronic correlations 
too, so that these systems are typical strongly correlated metals. Two types of correlated 
metals can be distinguished: (i) the doped Mott insulator and (ii) the bandwidth controlled 
correlated metal W ~ U. Both types will be considered in this paper. 

A state of the art tool to describe such correlated systems is the dynamical mean-field the- 
ory (DMFT) fr7i9i?i 1( ? . The characteristic features of correlated systems within the DMFT are 
the formation of incoherent structures, the so-called Hubbard bands, split by the Coulomb 
interaction U, and a quasiparticle (conduction) band near the Fermi level dynamically gen- 
erated by the local correlation o 6 i 7 i 8 i 9 i 10 . 
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Unfortunately, the DMFT is not useful to the study the "antiferromagnetic" scenario of 
pseudogap formation in strongly correlated metals. This is due to the basic approximation 
of the DMFT, which amounts to the complete neglect of non-local dynamical correlation 
effects. The aim of the present paper is to formulate an approach overcoming this difficulty 

The paper is organized as follows: In section |H] we present a derivation of the self- 
consistent generalization we call DMFT+Ek which includes short-ranged dynamical corre- 
lations to some extent. Section ITTT1 describes the construction of the k-dependent self-energy, 
and some computational details are presented in section IIV Al Results and a discussion are 
given in the sections HXl and El 



II. INTRODUCING LENGTH SCALE INTO DMFT: DMFT+S k APPROACH 

Basic shortcoming of traditional DMFT approac h 6 ' 7 ! 8 ! 9 ' 10 is the neglect of momentum 
dependence of electron self-energy. This approximation in principle allows for an exact 
solution of correlated electron systems fully preserving the local part of the dynamics in- 
troduced by electronic correlations. To include non-local effects, while remaining within 
the usual "impurity analogy", we propose the following procedure. To be definite, let us 
consider a standard one-band Hubbard model from now on. The extension to multi-orbital 
or multi-band models is straightforward. The major assumption of our approach is that the 
lattice and Matsubara "time" Fourier transform of the single-particle Green function can be 
written as: 

Gk(iw) = - — - s 1 r v ,. r , u = 7iT(2n + 1), (1) 

where E(ia>) is the local contribution to self-energy, surviving in the DMFT, while £k(iu;) is 
some momentum dependent part. We suppose that this last contribution is due to either elec- 
tron interactions with some "additional" collective modes or order parameter fluctuations, 
or may be due to non-local contributions within the Hubbard model itself. The assumed 
additive form of self-energy H(iu) + H^iiui) corresponds to neglect of possible interference 
of these local and non-local contributions. 

The self-consistency equations of our generalized DMFT+Sk approach are formulated 
as follows: 

1. Start with some initial guess of local self-energy E(ia>), e.g. E(ia>) = 0. 
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2. Construct H^iu) within some (approximate) scheme, taking into account interactions 
with collective modes or order parameter fluctuations which in general can depend on 
H{iuj) and \x. 

3. Calculate the local Green function 

= - ^ + ^ _ g(k) _ E( _ ^ ( . (2) 

4. Define the "Weiss field" 

^H^H + g^m. (3) 

5. Using some "impurity solver" to calculate the single-particle Green function for the 
effective Anderson impurity problem, defined by integral 

G d (r - T ') = ^- \ DcfDc iaCia (T)cUr') exp(-S eff ) (4) 

with effective action for a fixed site ("impurity") i 

Ses=- dn / dT 2 c icT (r 1 )g 1 (T 1 - r 2 )c+(r 2 ) + / drUn^i^ni^r) , (5) 
Jo Jo Jo 

Z cS = f Dcf a Dci a exp(— S e s), and (3 = T~ l . This step produces a new set of values 

G- d \iu). 

6. Define a new local self-energy 

sh^o^m-g^m. (6) 

7. Using this self-energy as "initial" one in step 1, continue the procedure until (and if) 
convergency is reached to obtain 

G u (iu) = G d {iu)- (7) 

Eventually, we get the desired Green function in the form of (Q), where Y,(iuj) and ^(iu) 
are those appearing at the end of our iteration procedure. A more detailed derivation of 
this scheme within a diagrammatic approach is given in the Appendix lAl 
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III. CONSTRUCTION OF K DEPENDENT SELF ENERGY 



For the momentum dependent part of the single-particle self-energy we concentrate on 
the effects of scattering of electrons from collective short-range SDW-like antiferromagnetic 
spin (or CDW-like charge) fluctuations. To calculate Sj c (za;) for an electron moving in 
the quenched random field of (static) Gaussian spin (or charge) fluctuations with dominant 
scattering momentum transfers from the vicinity of some characteristic vector Q, we use a 
slightly generalized version of the recursion procedure proposed in Refs . 4 ^ 11 and take into 
account all Feynman diagrams describing the scattering of electrons by this random field. 
Then the desired self-energy is given by 

£ k (iu;) = £ n=1 (iu;k) (8) 

with 

EJicuk) = A 2 ^ — - . (9) 

iuo + /i — zl{iLo) — e n \k) + inv n n — £ n+1 (iu;k) 

The quantity A characterizes the energy scale and k = £ _1 is the inverse correlation length 
of short range SDW (CDW) fluctuations, e n (k) = e(k + Q) and v n = I^+qI + l^k+ol f° r 
odd n while e n (k) = e(k) and v n = \v£\ + \v^\ for even n. The velocity projections and 
are determined by usual momentum derivatives of the "bare" electronic energy dispersion 
e(k). Finally, s(n) represents a combinatorial factor with 

s(n) = n (10) 

for the case of commensurate charge (CDW type) fluctuations with Q = (n/a, n/a)^. For 
incommensurate CDW fluctuations^ one finds 

{^rr- for odd n 
(11) 
| for even n. 

If we want to take into account the (Heisenberg) spin structure of interaction with spin 
fluctuations in "nearly antiferromagnetic Fermi-liquid" (spin-fermion (SF) model Ref.— ), the 
combinatorics of diagrams becomes more complicated. Spin-conserving scattering processes 
obey commensurate combinatorics, while spin-flip scattering is described by diagrams of 
incommensurate type ( "charged" random field in terms of Ref.—). In this model the recursion 
relation for the single-particle Green function is again given by 0, but the combinatorial 
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factor s(n) now acquires the following form 4 : 

( 2±2 for odd n 
s{n) = l 3 (12) 
I | for even n. 

Obviously, with this procedure we introduce an important length scale £ not present in 
standard DMFT. Physically this scale mimics the effect of short-range (SDW or CDW) cor- 
relations within fermionic "bath" surrounding the effective Anderson of the DMFT impurity. 
We expect that such a length-scale will lead to a competition between local and non-local 
physics. 

An important aspect of the theory is that both parameters A and £ can in principle 
be calculated from the microscopic model at hand. For example, using the two-particle 
selfconsistent approach of Ref-^ with the approximations introduced in Refs4»^, one can 
derive within the standard Hubbard model the following microscopic expression for A: 

A 2 = -U 2 k n iT > + < mi > -2 < n iT n a >1 = 



4 < X mi > 



U 2 < ( *T ~ ni l ) > 



where we consider only scattering from antiferromagnetic spin fluctuations. The different 
local quantities - spin fluctuation < S 2 >, density n and double occupancy < m-\ n ii > 
- can easily be calculated within the standard DMFT—. A detailed derivation of (fT^j) and 
computational results for A obtained by DMFT using quantum Monte-Carlo (QMC) to solve 
the effective impurity problem are presented in Appendix El A corresponding microscopic 
expressions for the correlation length £ = k," 1 can also be derived within the two-particle 
self-consistent approach"— However, we expect those results for £ to be less reliable, because 
this approach is valid only for relatively small (or medium) values of U/t. Thus, in the 
following we will consider both A and especially £ as some phenomenological parameters to 
be determined from experiments. 
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IV. RESULTS AND DISCUSSION 



A. Computation details 

In the following, we want to discuss results for a standard one-band Hubbard model on 
a square lattice. With nearest (t) and next nearest (t 1 ) neighbour hopping integrals the 
dispersion then reads 

e(k) = — 2t(cos k x a + cos k y a) — 4t' cos k x a cos k y a , (14) 

where a is the lattice constant. The correlations are introduced by a repulsive local two- 
particle interaction U. We choose as energy scale the nearest neighbour hopping integral t 
and as length scale the lattice constant a. 

For a square lattice the bare bandwidth is W = 8t. To study a strongly correlated 
metallic state obtained as doped Mott insulator we use U = 40t as value for the Coulomb 
interaction and a filling n = 0.8 (hole doping). The correlated metal in the case of W > U 
is realized via U = 4t and two fillings: half-filling (n = 1.0) and n = 0.8 (hole doping). As 
typical values for A we choose A = t and A = 2t (actually as approximate limiting values 
- cf. Appendix |B} and for the correlation length £ = 2a and £ = 10a (motivated mainly 
by experimental data for cuprates^). 

The DMFT maps the lattice problem onto an effective, self-consistent impurity defined 
by Eqs. (@J)-(JH1). In our work we employ as "impurity solvers" two reliable numerically 
exact methods — quantum Monte-Carlo (QMC) 13 and numerical renormalization group 
(NRG) 15 i 16 . Calculations were done for the case t' = and t'/t=-0A (more or less typical 
for cuprates) at two different temperatures T = 0.088t and T = 0.356i (for NRG compu- 
tations)^ . QMC computations of double occupancies as functions of filling were done at 
temperatures T = O.lt and T = 0.4t 22 . 

B. Generalized DMFT+Sk approach: densities of states 

Let us start the discussion of our results obtained within our generalized DMFT+Sk 
approach with the densities of states (DOSs) for the case of small (relative to bandwidth) 
Coulomb interaction U = 4i with and without pseudogap fluctuations. As already discussed 
in the introduction, the characteristic feature of the strongly correlated metallic state is the 

7 



coexistence of lower and upper Hubbard bands split by the value of U with a quasiparticle 
peak at the Fermi level. Since at half-filling the bare DOS of the square lattice has a 
Van-Hove singularity at the Fermi level (t' = 0) or close to it (in case of t'/t = —0.4) one 
cannot treat a peak on the Fermi level simply as a quasiparticle peak. In fact, there are two 
contributions to this peak from (i) the quasiparticle peak appearing in strongly correlated 
metals due to many-body effects and (ii) the smoothed Van-Hove singularity from the bare 
DOS 17 . In Figs, ^and |21 we show the corresponding DMFT(NRG) DOSs without pseudogap 
fluctuations as black lines for n = 1 and n = 0.8 for both bare dispersions t'/t = —0.4 
(left panels) and for t' = (right panels) for two different temperatures T = 0.356t (middle 
panels) and T = 0.088t (upper and lower panels). The remaining curves in Figs. Q and 121 
represent results for the DOSs with non-local fluctuations switched on. The fluctuation 
amplitudes are taken as A = It (upper and and middle panels) and A = t (low panels). For 
all sets of parameters one can see that the introduction of non-local fluctuations into the 
calculation leads to the formation of pseudogap in the quasiparticle peak. 

The behaviour of the pseudogaps in the DOSs has some common features. For example, 
for t'=0 at half-filling (Fig. ^ right column) we find that the pseudogap is most pronounced. 
For n = 0.8 (Fig. |21 right column) the picture is almost the same but slightly asymmetric. 
The width of the pseudogap (the distance between peaks closest to Fermi level) appears to 
be of the order of ~ 2A here. Decreasing the value of A from 2t to t leads to a pseugogap 
that is correspondingly twice smaller and in addition more shallow. When one uses the 
combinatorial factors corresponding to the spin-fermion model (Eq. (jl2j) ). we find that the 
pseudogap becomes more pronounced than in the case of commensurate charge fluctuations 
(combinatorial factors of Eq. (JUJ). The influence of the correlation length £ can be seen is 
also as expected. Changing £ _1 from £ _1 = 0.1 to £™ 1 = 0.5, i.e. decreasing the range of the 
non-local fluctuations, slightly washes out the pseudogap. Also, increasing the temperature 
from T = 0.088t to T = 0.356t leads to a general broadening of the structures in the 
DOSs. These observations remain at least qualitatively valid for t'/t = —0.4 (Figs. Hand 
left columns) with an additional asymmetry due to the next-nearest neighbour hopping. 
Noteworthy is however the fact that for t'/t = —0.4 and £ _1 = 0.5 the pseudogap has 
almost disappeared for the temperatures studied here. An also very remarkable point is the 
similarity of the results obtained with the generalized DMFT+Ek approach with U — At 
(smaller than the bandwidth W) to those obtained earlier without Hubbard-like Coulomb 
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interactions*^. 

Let us now consider the case of a doped Mott insulator. The model parameters are again 
t' = and t'/t = —0.4 with filling n = 0.8, but the Coulomb interaction strength is set 
to U = AOt for our DMFT+Sk calculations. Characteristic features of the DOS for such a 
strongly correlated metal are a strong separation of lower and upper Hubbard bands and 
a Fermi level crossing by the lower Hubbard band (for non-half-filled case). Without non- 
local fluctuations the quasi-particle peak is again formed at the Fermi level; but now the 
upper Hubbard band is far to the right and does not touch the quasiparticle peak (as it was 
for the case of small Coulomb interactions). DOSs without non-local fluctuations are again 
presented as black lines on Fig. HJ (£' = 0) and Fig. |3] (t'/t = —0.4). 

With rather strong non-local fluctuations A = 2t, a pseudogap appears in the middle of 
quasiparticle peak. In addition we observe that the lower Hubbard band is slightly broadened 
by fluctuation effects. Qualitative behaviour of the pseudogap anomalies is again similar to 
those described above for the case of U = At, e.g. a decrease of £ makes the pseudogap less 
pronounced, reducing A from A = 2t to A = t narrows of the pseudogap and also makes 
it more shallow etc.. Note that for the doped Mott-insulator we find that the pseudogap is 
remarkably more pronounced for the SDW-like fluctuations than for CDW-like fluctuations. 

There are, however, quite clear differences to the case with U = At. For example, the 
width of the pseudogap appears to be much smaller than 2A, which we attribute to the fact 
that the quasiparticle peak itself is actually rather narrow now. 



C. Generalized DMFT+Sk approach: spectral functions A(uj, k) 

In the previous subsections we discussed the densities of states obtained self-consistently 
by the DMFT+Sk approach. Once we get a self-consistent solution of the DMFT+Ek 
equations with non-local fluctuations we can of course also compute the spectral functions 
A(u,k) 

1 1 

A(uj,k) = — Im — ■ — — —— — ^r-r^> ( 15 ) 

where self-energy £(u;) and chemical potential // are calculated self-consistently as described 
in Sec. |H] To plot A(u,\t) we choose k-points along the "bare" Fermi surfaces for different 
types of lattice spectra and fillings. In Fig. El one can see corresponding shapes of these 
"bare" Fermi surfaces (presented are only 1/8-th parts of the Fermi surfaces within the first 
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quadrant of the Brillouin zone). 

In the following we concentrate mainly on the case U = At and filling n = 0.8 (Fermi 
surface of Fig. E(a)). The corresponding spectral functions A(u,k) are depicted in Fig. [7| 
When if ft = —0.4 (upper row), the spectral function close to the Brillouin zone diagonal 
(point B) has the typical Fermi-liquid behaviour, consisting of a rather sharp peak close to 
the Fermi level. In the case of SDW-like fluctuations this peak is shifted down in energy 
by about — 0.5t (left upper corner). In the vicinity of the "hot-spot" the shape of A(u, k) 
is completely modified. Now A(u, k) becomes double-peaked and non-Fermi-liquid-like. 
Directly at the "hot-spot", A(u;,k) for SDW-like fluctuations has two equally intensive 
peaks situated symmetrically around the Fermi level and split from each other by ~ 1.5 A 
Refs.^. For commensurate CDW-like fluctuations the spectral function in the "hot-spot" 
region has one broad peak centred at the Fermi level with width ~ A. Such a merging 
of the two peaks at the "hot-spot" for commensurate fluctuations was previously observed 
in RefA However close to point A this type of fluctuations also produces a double-peak 
structure in the spectral function. 

In Fig. El we show spectral functions for the case of U = At at half-filling (n = 1) 
(Fermi surface of Fig. El (c), (d)). For tl jt = —0.4 (upper row of Fig. EJ) these are similar 
to those just discussed for n = 0.8. However, the pseugogap is more pronounced here and 
remains open even close to the antiferromagnetic Fermi surface boundary (point B) for SDW 
fluctuations. 

In the lower panel of Fig. Owe show spectral functions for 20% hole doping (n = 0.8) and 
the case of t! — (Fermi surface from Fig. Efb)). Since the Fermi surface now is close to 
the antiferromagnetic zone boundary, the pseudogap anomalies are rather strong and almost 
non-dispersive along the Fermi surface. At half-filling (Fig. [?1 lower panel) and for if — 
the Fermi surface (Fig. 03(d)) actually coincides with antiferromagnetic zone boundary. In 
this case the the whole Fermi surface is the "hot-region" (perfect "nesting"). The spectral 
functions are now symmetric around the Fermi level. For SDW-like fluctuations there are 
two peaks split by ~ 1.5 A. Again, CDW-like fluctuations give just one peak centred at the 
Fermi level with width ~ A. 

For the case of a doped Mott insulator (U = AOt, n = 0.8), the spectral functions obtained 
by the DMFT+S k approach are presented in Fig. |H1 Qualitatively, the shapes of these 
spectral functions are similar to those shown on Fig. As was pointed out above, the strong 
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Coulomb correlations lead to a narrowing of the quasiparticle peak and a corresponding 
decrease of the pseudogap width. One should also note that in contrast to U — At the 
spectral functions are now less intensive, because part of the spectral weight is transferred 
to the upper Hubbard band located at about 401 

Using another quite common choice of k-points we can compute A(u, k) along high- 
symmetry directions in the first Brillouin zone: r(0, 0) — X(7r, 0) — M(jr, tt) — r(0, 0). The 
spectral functions for these k-points are colelcted in Figs. El and ^2 for the case of 
SDW-like fluctuations. For all sets of parameters one can see a characteristic double peak 
pseudogap structure at the X point and in the middle of M — T direction. A change of the 
filling leads mainly to a rigid shift of spectral functions with respect to the Fermi level. 

With the spectral functions we are now of course in a position to calculate angle resolved 
photoemission spectra (ARPES), which is the most direct experimental way to observe 
pseudogap in real compounds. For that purpose, we only need to multiply our results for 
the spectral functions with the Fermi function at temperature T = 0.088t. the resulting 
DMFT+Ek ARPES spectra are presented in Figs. fT2| II HI and ITU One should note that for 
t'/t = —0.4 (upper panel of Figs. H2J CHI and Ej) as k goes from point "A" to point "B" 
the peak situated slightly below the Fermi level changes its position and moves down in 
energy. Simultaneously it becomes more broad and less intensive. The dotted line guides 
the motion of the peak maximum. Also at the "hot-spot" and further to point B one can 
see some signs of the double-peak structure. Such behaviour of the peak in the ARPES is 
rather reminiscent of those observed experimentally in underdoped cuprates^*^. 

V. CONCLUSION 

To summarize, we propose a generalized DMFT+Ek approach, which is meant to take 
into account the important effects due to non-local correlations in a systematic, but to 
some extent phenomenological fashion. The main idea of this extension is to stay within 
a usual effective Anderson impurity analogy, and introducing length scales due to non- 
local correlation via the effective medium ("bath") appearing in the standard DMFT. This 
becomes possible by incorporating scattering processes of fermions in the "bath" from non- 
local collective SDW-like antiferromagnetic spin (or CDW-like charge) fluctuations. Such 
a generalization of the DMFT allows one to overcome the well-known shortcoming of k- 
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independence of self-energy of standard DMFT. It in turn opens the possibility to access 
the physics of low-dimensional strongly correlated systems, where different types of spatial 
fluctuations (e.g. of some order parameter), become important in a non-perturbative way at 
least with respect to the important local dynamical correlations. However, we must stress 
that our procedure in no way introduces any kind of systematic l/<i-expansion, being only 
a qualitative method to include length scale into DMFT. 

In our present study we addressed the problem of pseudogap formation in the strongly 
correlated metallic state. We showed evidence that the pseudogap appears at the Fermi 
level within the quasiparticle peak, introducing a new small energy scale of the order of A 
in the DOSs and spectral functions A(u, k). 

Let us stress, that our generalization of DMFT leads to non-trivial and in our opinion 
physically sensible k-dependence of spectral functions. In contrast, in a recent work by 
Haule et a/.— the extended DMFT approach was used to demonstrate pseudogap formation 
in DOS due to dynamic Coulomb correlations only. However, within the approach of Ref.— 
there is no way to obtain k-dependence of spectral functions beyond that originating from 
the bare electronic energy dispersion which is actually observed in experiments. Of course, 
similar results and observations were in recent years also made using the cluster mean-field 
theories^. The major advantage of our approach over these cluster mean-field theories 
is, that we stay in an effective single- impurity picture. This means that our approach is 
computationally much less expensive and therefore also easily generalizable to multi-orbital 
systems. 
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APPENDIX A: DERIVATION OF GENERALIZED DMFT+S k APPROACH 



In this appendix we present a derivation of the generalized DMFT+Ek scheme for the 
Hubbard model 

H = - *ii4r c *r + U n i\ n ih ( Al ) 
ij,cr i 

using a diagrammatic approach. The single-particle Green function in Matsubara represen- 
tation is as usual given by 

G ^ = -T ( L w (A2) 

ilu + fi — e(k) — h(iu, k) 

To establish the standard DMFT one invokes the limit of infinite dimensions d — > oo. In 
this limit only local contributions to electron self-energy survive^, i.e. E^- — > 5yEjj or, in 
reciprocal space, E(iu;,k) — > H(iu). 

In Fig. E3we show examples of "skeleton" diagrams for the local self - energy, contribut- 
ing in the limit of d — > oo. The complete series of these and similar diagrams defines the 
local self - energy as a functional of the local Green function 

E = F[Gn] , (A3) 

where 

G*M = — * w . r - (A4) 

N ilj + fi — e{K) — H[zu) 

One then defines the "Weiss field" 

^M^M + G^M (A5) 

which is used to set up the effective impurity problem with an effective action given by 
(JHJ). Via Dyson's equation the Green function (jlj) for this effective impurity problem can be 
written as 



G dM = ^7—. (A6) 



1 

and the "skeleton" diagrams for self-energy E^ are just the same as shown in Fig. with 
the replacement Gu — * Gd- Thus we get 



S d = F[G d ], (A7) 



where F is the same functional as in (|A3|) . The two equations (jA6|) and (|A7|) define both 
and E d for a given "Weiss field" Q . On the other hand, for the local E and Gu of the 
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initial (Hubbard) problem we have precisely the same pair of equations, viz (jA3|) and (|A5|) . 
and Go i n both problems is just the same, so that 



Gu — Gd- 



(A8) 



Thus, the task of finding the local self-energy of the (d — > oo) Hubbard model is eventually 
reduced to the calculation of the self-energy of an effective quantum impurity problem 
defined by effective action of Eq. (JHJ). 

Consider now non - local contribution to the self - energy. If we neglect interference 
between local and non-local contributions (as given e.g. by the diagram shown in Fig lToT bD. 
the full self-energy is approximately determined by the sum of these two contributions. 
"Skeleton" diagrams for the non-local part of the self-energy, Sk(^), are then those shown 
in Fig. IToT a). where the full line denotes the Green function of Eq. while broken 
lines denote the interaction with static Gaussian spin (charge) fluctuations. 

The local contribution to the self-energy is again defined by the functional (|A3J) via the 
local Green function Ga, which is now given by (J2J). Introducing again a "Weiss field" via 
(IA5|) and repeating all previous arguments, we again reduce the task of finding the local 
part of the self-energy to the solution of an "impurity" problem with an effective action (j3J) . 

To determine the non-local contribution H^iu) we first introduce 



as the "bare" Green function for electron scattering by static Gaussian spin (charge) fluctua- 
tions. The assumed static nature of these fluctuations allows to use the method of Refs . 4l5lU 
and the calculation of the non-local part of the self-energy H^iu) reduces to the recursion 
procedure defined by Eqs. (jHI) and Q- The choice of the "bare" Green function Eq. (|A9|) 
guarantees that the Green function "dressed" by fluctuations G^^uj) = Qq^^) — £k(^), 
which enters into the "skeleton" diagrams for T>k(iu), just coincides with the full Green 
functions G^iu). 

Thus we obtain a fully self-consistent scheme to calculate both local (due to strong single- 
site correlations) and non-local (due to short-range fluctuations) contributions to electron 
self-energy. 



1 1 



(A9) 



G^^lu) + Ek(iu) iuj + fi — e(k) — luiiu) 
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APPENDIX B: A IN THE HUBBARD MODEL. 



In this Appendix we derive the explicit microscopic expression for pseudogap amplitude 
A given in (|13|). Within the two-particle self-consistent approach of Ref.— , valid for medium 
values of U, and neglecting charge fluctuations, we can write down an expression for the 
electron self-energy of the form used in (0), with 

E a {iu) = Un_ a (Bl) 

as the lowest order local contribution due to the on-site Hubbard interaction, surviving in the 
limit of d —>■ oo, and exactly accounted for in DMFT (with all higher-order contributions). 
Non-local contribution to the self-energy (vanishing for d — > oo) due to interaction with 
spin-fluctuations then lead to the expression 

U T . 

Ek(«w) = -TT7 UwXwfa* u rn)G (k + q, iu + iv m ) , (B2) 



4 N 

m q 



where 



U v = g n (0)U, ffTl(O) = ~ < n Tj l > Z ( B3 ) 

< n iT >< nn > 

with < r? a >=< n a > and < >=< >= \n in the paramagnetic phase. For the 
dynamic spin susceptibility Xsp{<l, v m) we use the standard Ornstein-Zernike forrn^, similar 
to that used in spin-fermion model^, which describes enhanced scattering with momenta 
transfer close to antiferromagnetic vector Q = (n/a, n/a). With these approximations, we 
can write down the following expression for the non-local contribution to the self-energy^: 

1 



1 T x 



N ^— ' ^— ' ' zu; + zz/ m + /i — e(k + q) 



q 



1 



™ ^ „ zu; + /i — e(k + q) 

1 



q iu + fi- e(k + q) 



A 



2 



iw + /i - e(p + Q) + i(|t£+ q | + | ^ +Q |) resigns ' 
Here we have introduced the static form factor 5 



(B4) 



s{q) (& - 1) 2 + r 2 - q!) 2 + e~ 2 (B5) 
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and the squared pseudogap amplitude 

1 T . 

A 2 = -^UU sp — Yl Xs p^ Vm ) = 

m q 

= hjU sp [< n iT > + < > -2 < n^n^ >] = 

= ~UU sp ~ < S 2 >, (B6) 

where we have used the exact sum-rule for the susceptibility 4,12 . Taking into account ()B3j) 
we immediately obtain (fT3|). 

Actually, the approximations made in (|B4|) and ()B5|) allow for an exact summation of the 
whole Feynman series for electron interaction with spin-fluctuations, replaced by the static 
Gaussian random field. Thus generalizing the one-loop approximation (|B4j) eventually leads 
to the basic recursion procedure given in (JHJ), © Refs4£. 

Using the DMFT(QMC) approach we computed occupancies < >, < > and double 
occupancies < n^Tty > required to calculate the pseudogap amplitude A of Eq. (jB6|) In Fig. 
El the corresponding values of A are presented. One can see that A grows when the filling 
goes to n = 1. While U approaches 8t (the value of the bandwidth for a square lattice) A as 
a function of n grows monotonically. When U becomes larger than W = 8t (when a metal- 
insulator transition occurs) one can see a local minimum for n = 0.9, which becomes more 
pronounced with further increase of U. For t'/t = —0.4 and both temperatures the scatter 
of A values is smaller than for the case of t' — 0. Also A has a rather weak temperature 
dependence. All values of A lie in the interval ~ 0.75£-i-2£. Therefore, for our computations 
we took only two characteristic values of A = t and A = 2t. 
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FIG. 1: Comparison of DOSs obtained from DMFT(NRG)+Sj c calculations for different combi- 
natorical factors (SF — spin-fermion model, commensurate), inverse correlation lengths in 
units of the lattice constant, temperatures (T) and values of the pseudogap potential (A). Left 
column corresponds to t'/t = —0.4, right column to t' = 0. In all graphs the Coulomb interaction 
is U = At and n = 1. The Fermi level corresponds to zero. 
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FIG. 2: Comparison of DOSs obtained from DMFT(NRG)+£k calculations for a filling n = 0.8, 
other parameters as in Fig. ^ 
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FIG. 3: Comparison of DOSs obtained from DMFT(NRG)+S k calculations for t'/t = -0.4, T = 
0.088t, U = 40t and filling n = 0.8. 





FIG. 4: Comparison of DOSs obtained from DMFT(NRG)+S k calculations for if = 0, other 
parameters as in Fig. |31 
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FIG. 5: 1/8-th of the bare Fermi surfaces for the different occupancies n and combinations (t,t') 
used for the calculation of spectral functions A(k,uj). The diagonal line corresponds to the anti- 
ferromagnetic Brillouin zone boundary at half-filling for a square lattice with nearest-neighbours 
hopping only. The full circle marks the so-called "hot-spot" . 
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FIG. 6: Spectral functions A(k,uj) obtained from DMFT(NRG)+Sk calculations along the direc- 
tions shown in Fig. [5J Model parameters were chosen as U = 4i, n = 1.0, A = 2t, = 0.1 and 
temperature T = 0.088t. The "hot-spot" k-point is marked as fat dashed line. The Fermi level 
corresponds to zero. 
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FIG. 7: Spectral functions A(k, uj) obtained from the DMFT(NRG)+$]k calculations for n = 0.8, 
other parameters as in Fig. HO 



_ B 
3 

M 
< 

a 

■S Al 
o 
a 

a 

00 



A 



r i 1 i 1 i 1 i 1 i 

Commensurate 
t'/t=-0.4 




For all graphs 
A=2t, U=40t, T=0.088t, 11=0.8, if '=0.1 




-3-2-1 1 2 3 4 -2-1 1 2 3 4 5 

Energy in units of t 



FIG. 8: Spectral functions A(k, u) obtained from the DMFT(NRG)+I]k calculations for U 
other parameters as in Fig. [7| 




FIG. 9: Spectral functions A(k, uj) obtained from the DMFT(NRG)+S k calculations along high- 
symmetry directions of first Brilloin zone r(0, 0)— X(7r, 0)— M(ir, tt)— T(0, 0), SF combinatorics (left 
row) and commensurate combinatorics (right column). Other parameters are U = 4t, n = 1.0, 
A = 2t, ^ 1 = 0.1 and temperature T = 0.088t. The Fermi level corresponds to zero. 
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FIG. 10: Spectral functions A(k, uS) obtained from the DMFT(NRG)+Sk calculations along high- 
symmetry lines in the 2D Brillouin zone for n = 0.8, other parameters as in Fig. |SJ 
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FIG. 11: Spectral functions A(k,uj) obtained from the DMFT(NRG)+Sk calculations along high- 
symmetry lines in the 2D Brillouin zone for U = 40t, other parameters as in Fig. I1UI 




FIG. 12: ARPES obtained from the DMFT (NRG) + S k calculations for U = At and n = 1.0 along 
the lines in the first BZ as depicted by Fig. |SJ all other parameters as in Fig. 
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FIG. 13: ARPES obtained from the DMFT (NRG) + S k calculations for U = At and n = 0.8 along 
the lines in the first BZ as depicted by Fig. |SJ all other parameters as in Fig. [7| 
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FIG. 14: ARPES obtained from the DMFT (NRG) + S k calculations for U = 40t and n = 0.8 along 
the lines in the first BZ as depicted by Fig. |SJ all other parameters as in Fig. |HJ 
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FIG. 15: Local "skeleton" diagrams for the DMFT self-energy S. Wavy lines represent the local 
(Hubbard) Coulomb interaction U, full lines denote the local Green function Ga. 
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FIG. 16: Typical "skeleton" diagrams for the self-energy in the DMFT+Sk approach. The first 
two terms are DMFT self-energy diagrams; the middle two diagrams show contributions to the 
non-local part of the self-energy from spin fluctuations (see sectior JTTl|) represented as dashed lines; 
the last diagram (b) is an example of neglected diagramms leading to interference between the 
local and non-local parts. 




FIG. 17: Filling dependence of the pseudo gap potential A calculated with DMFT(QMC) for 
varying Coulomb interaction (U) and temperature (T) on a two-dimensional square lattice with 
two sets of (t, t'). 



